Fisher information under decoherence in Bloch representation 
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The dynamics of two variants of quantum Fisher information under decoherence are investigated 
from a geometrical point of view. We first derive the exphcit formulas of these two quantities for 
a single qubit in terms of the Bloch vector. Moreover, we obtain analytical results for them under 
three different decoherence channels, which are expressed as affine transformation matrices. Using 
the hierarchy equation method, we numerically study the dynamics of both the two information in a 
dissipative model and compare the numerical results with the analytical ones obtained by applying 
the rotating- wave approximation. We further express the two information quantities in terms of the 
Bloch vector for a qudit, by expanding the density matrix and Hermitian operators in a common set 
of generators of the Lie algebra 5u{d). By calculating the dynamical quantum Fisher information, 
we find that the collisional dephasing significantly diminishes the precision of phase parameter with 
the Ramsey interferometry. 

PACS numbers: 03.65.Yz, 03.65.Ta, 03.67.-a 



I. INTRODUCTION 

Quantum Fisher information (QFI), which is one of 
the most important quantities for both quantum estima- 
tion theory and quantum information theory, has been 
widely studied [1-9]. In the field of quantum estimation, 
the main task is to determine the value of an unknown 
parameter labeling the quantum system, and a primary 
goal is to enhance the precision of the resolution [10-20]. 
The inverse of the QFI provides the lower bound of error 
of the estimation [21, 22]. Hence, how to increase the QFI 
become the key problem to be solved. Moreover, the QFI 
can be used to measure the statistical distinguishability 
on the space of the density operators in the quantum in- 
formation geometry [23, 24]. Recently, the QFI flow was 
proposed as a quantitative measure of the information 
flow and provides a novel perspective on observing the 
non-Markovian behavior in open quantum systems [25]. 

There are several variants of quantum versions of 
Fisher information, among which the one based on the 
symmetric logarithmic derivative (SLD) operator has 
been used most widely and possesses many good proper- 
ties [21, 22], such as convexity, remaining invariant under 
the unitary evolution, and the total amount of the QFI 
equivalent to the summation of the QFIs of all subsystem 
being uncorrelated. In this paper, we also study another 
variant of the quantum version of the Fisher informa- 
tion, which is closely related with the skew information 
[26] . The skew information was proposed to measure the 
amount of information that a quantum state contains 
with respect to the observable which does not commute 
with additive conserved quantities, such as Hamiltoni- 
ans, momenta [26]. It is also used to measure the quan- 
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tum uncertainty [27] and quantify quantum correlations 
in bipartite state in a recent work [28]. In the following, 
we shall show that the two QFIs formally possess many 
similar features [1], but they are different not only for 
the concrete expressions of the definitions, but also for 
their applications. As is mentioned previously, the QFI 
based on SLD operator is mainly applied in the quan- 
tum metrology, however, the variant version of the QFI 
plays an important role in quantum state discrimination 
[29, 30]. 

From the information theory, both two different QFIs 
characterize the content of information contained in the 
quantum systems. A crucial property of these two infor- 
mation quantities is that they decrease monotonically un- 
der the completely positive and trace-preserving (CPT) 
maps [1, 31-33]. The monotonicity property manifests 
the information loss under the CPT map. Just as the in- 
troduction of the QFI flow [25], can we observe the non- 
Markovian properties from the other information quan- 
tities perspective? 

Here, we will address this problem by calculating these 
two QFIs for a single qubit and derive the explicit formu- 
las in the Bloch representation, which greatly facilitates 
the computing of these two quantities. The main results 
in this paper are that the dynamical QFIs, in the pres- 
ence of decoherences, are analytically solved. Here, the 
irreversible processes are modeled via three decoherence 
channels [7, 34-36]: phase-damping channel (PDC), de- 
polarizing channel (PDC), and the generalized amplitude 
damping channel (GADC) [37]. The analytical results 
for the two information quantities under those channels 
are obtained. We will show that the values of the two 
QFIs monotonically decrease with time, apart from the 
isolated case that the QFI based on the SLD operator 
about the amplitude parameter 6 remain invariant un- 
der the PDC. In order to further identify the behaviors 
of the two QFIs subject to quantum noise, we discuss a 
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simple model of a two-level system coupled to a reservoir 
with a Lorentzian spectral density [38, 39]. By using the 
hierarchy equation method [40-47], we numerically ana- 
lyze the dynamics of the two QFIs, and compare these 
with the analytical solutions by using the rotating-wave 
approximation (RWA). We further generalize the results 
to the qudit system. Meanwhile, we verify that they are 
also applicable for the 7V-qubit system with symmetry 
exchange. For the sake of clarity, we calculate the dy- 
namical quantum Fisher information in the presence of 
collisional dephasing. 

This paper is organized as follows. In Sec. II, we first 
review two different definitions of the QFI, and give the 
explicit formulae for the QFIs for a single qubit sys- 
tem. In Sec. Ill, we obtain the analytical results for the 
two quantities under three different decoherence chan- 
nels, and the numerical results are given. Moreover, in 
Sec. IV(a), we generalize the expressions of the two QFIs 
for a qudit system, and the QFI in a noisy environment 
for an 7V-qubit system is discussed in Sec. IV(b). Finally, 
the conclusion are given in Sec. V. 



II. FISHER INFORMATION 

In this section, we briefly summarize two variants of 
definition of the QFIs [26, 48], which are referred to as 
two different extensions from the classical Fisher infor- 
mation. We also discuss the relations between the QFI 
and the Bures distance [24, 49-51] as well as the QFI and 
the Hellinger distance [52]. In the Bloch representation, 
we derive the explicit formulas of these two information 
quantities for the single qubit system. 



A. Fisher information and Bures distance 

The classical Fisher information, originating from the 
statistical inference, is a way of measuring the amount 
of information that an observable random variable X 
carries about an unknown parameter A. Suppose that 
{pi (A) , A G M}^^ is the probability density conditioned 
on the fixed value of the parameter A = A* with measure- 
ment outcomes {xi]. The classical Fisher information is 
defined as 



Fx = Y,Pi (A) 



a In Pi (A) 



(1) 



which characterizes the inverse variance of the asymp- 
totic normality of a maximum-likelihood estimator. Here 
we have assumed that the observable X is a discrete vari- 
able. If it is continuous, the summation in Eq. (1) should 
be replaced by an integral. 

The quantum analog of the Fisher information is for- 
mally generalized from Eq. (1) and defined as 



in terms of the symmetric logarithmic derivative (SLD) 
operator L^, which is a Hermitian operator determined 

by 



dxP\ = \{P\,Lx). 



(3) 



where d\ = ^ and { • , • } denotes the anticommu- 
tator. By diagonalizing the density matrix as p\ = 
Y.iQi\^i) {"^^^ associated with > and Y^iQi = 
the elements of the SLD operator are completely defined 
under the condition qi + Qj ^ 0. Therefore, Eq. (2) can 
be expressed as 



md,^,)\\ (4) 



where the first and the second summations involve sums 
over all qi' ^ and Qi -\- Qj 7^ 0, respectively. In Eq. (4), 
the first term is equal to the classical Fisher information 
of Eq. (2), which is called the classical term, and the 
second term is called the quantum term. For pure states, 
Eq. (4) reduces to 



(5) 



An essential feature of the QFI is that we can obtain the 
achievable lower bound of the mean-square error of the 
unbiased estimator for the parameter A, i.e., the so-called 
quantum Cramer-Rao (QCR) theorem: 



Var(A) > 



1 



(6) 



where Var (•) denotes the variance, A denotes the unbi- 
ased estimator, and u represents the number of repeated 
experiments. 

As shown in a seminal work, see Ref. [24], they ob- 
served that the estimability of a set of parameters param- 
eterizing the family of the quantum states {px}^ which is 
characterized by the QFI, is naturally related to the dis- 
tinguishability of the states on the manifold of the quan- 
tum states, which is measured by the Bures distance. 
They also proved that the QFI is simply proportional to 
the Bures distance [24] 

Dl[p{X),p{X + dX)] = ^TxdX\ (7) 

The Bures distance measures the distance between two 
quantum states, and is defined as [49-51] 



Di {p,a) := 2 (1 - TrV p'/'ap^/') , 



(8) 



(2) 



where the second term in the bracket is the so-called 
Uhlmann fidelity [50]. Meanwhile, the explicit formula 
of the QFI for the two-dimensional density matrices is 
obtained by [53] 

J^x = Tr{dxpf + ^Tv{pdxpf. (9) 
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In the Bloch sphere representation, any qubit state can 
be written as 



(10) 



where u? = {ujx^ cj^, cj^)^ is the real Bloch vector and a = 
{^xi ^z) denotes the Pauli matrices. Apparently, the 
eigenvalues of the density operator are (1 zb(x;)/2, with 
the length of the Bloch vector u = \u)\. Here, the Blcoh 
vector satisfies < 1, the equality holds for pure states. 
In the Bloch representation, can be represented as 
follows 



w = 1. 



(11) 



The first line of the above equation is only applicable for 
the mixed states, which can be straightly obtained by 
substituting Eq. (10) into Eq. (9). 

For pure states, we have equation = p. Taking 
differential on both side of this equation with respect to 
A, one gets 

dxp = dxp" = p (dxp) + (dxp) p. (12) 

The SLD operator is given as 

L = 2dxp, (13) 

by comparing Eq. (12) with Eq. (3). Substituting 
Eqs. (10) and (13) into Eq. (2), and using the relation 



Tr[(a-6-) {b '&)]=2a b 



(14) 



finally yield the second line of Eq. (11), i.e., Tx for pure 
states is the norm of the derivative of the Bloch vector. 



B. Fisher information and Hellinger distance 

By extending Eq. (1) into the quantum regime in a 
diff"erent way, we will obtain a variant QFI [54]. Eq. (1) 
can be equivalently expressed as 



F,=4^(aAVMA))'. 



(15) 



By straightly replacing the summation by a trace, the 
probability pi (A) by a density matrix px, and the differ- 



G,- 



with G 



ential dx by the inner differential dx' = i 

being a fixed self-adjoint operator and [ •, • ] denoting the 
commutator. One obtains the following equation 



lx:=4Tr (Sav^J =-4Tr 



In the particular case, suppose that px = 
i.e., Px satisfies the Landau- von Neumann equation 



(16) 



pc- , 



'^dxpx = G^px ' Note that here we define Eq. (16) as 

a paradigmatic version of the QFI. Actually, Eq. (16) is 
the so-called skew information 

XwY 



Px 



(17) 



which is introduced by Wigner and Yanase [26], with 
ignorance of a negligible constant number here, i.e. X = 
8XwY- The skew information (SI) is a measure of the 
information contained in a quantum state px with respect 
of a fixed conserved observable G. 

By inserting the spectrum decomposition p = 
Qi into Eq. (16), one obtain 



{d\Qi') 



Wj) \{'il'i\dxi}j)Y 



(18) 

where the first summation is over all 7^ 0, the same 
requirement as in Eq. (18). Comparing with Eq. (4), the 
classical terms in Eqs. (4) and (18) are the same, but the 
quantum terms are different. For pure states, Eq. (18) 
reduces to 



{dxMdx^)-mdxm 



(19) 



which is twice as much as the QFI of Eq. (2) given in 
Eq. (5). The factor-of-2 difference between Fx and Xx 
result from their different coefficients (or weights) of the 
quantum terms in Eqs. (4) and (18). For pure states, 
the classical terms vanish, and only the quantum terms 
contribute. Then we obtain Eqs. (5) and (19). 

Similar to the relation between the QFI of Eq. (2) and 
the Bures distance given in Eq. (7), there formally ex- 
ists the same relation between the QFI of Eq. (17) and 
quantum Hellinger distance [52]: 



D2 [p{X),p{X^d\)] = -Ixd\^ 



(20) 



The latter is the quantum version of the classical 
Hellinger distance, which measures the distance between 
two probability distributions. In the space of the quan- 
tum state, the quantum Hellinger distance is defined as 



^QH {P,^) := 2 (1 - Trv^v^) 



(21) 



where the last term in the bracket is called quantum aflin- 
ity [52]. 

For a 2 X 2 density matrix, we derive the explicit for- 
mula of the QFI Xx as follows: 

Xx = aTiidxpf -f3TT{pdxp)\ (22) 
where the coefficients are determined by 



1 



1 -4detp 
1 



4(l-2detp) 

(rT27asfpJ 



1 - 4 det p V 1 + 2 Vdet p det p J ' 



(23) 
(24) 
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In the Bloch representation, Xx can be represented as 
follows 

X, JtOT+®"<-*")'- (25) 

\2\dM\ ^ = h 

where the coefficient is given by 




The first line of Eq. (25) is only applicable for mixed 
states, which can be directly derived by substituting 
Eq. (10) into Eqs. (22), (23) and (24). 
For pure states, we have 

^p = p. (27) 

With Eqs. (14) and (27), one can obtain the second line 
of Eq. (25) from Eq. (16), i.e., Xx for pure states is the 
norm of the derivative of the Bloch vector up to a factor 
of 2. 



C. Example 

Having obtained Eqs. (11) and (25), we consider a sim- 
ple example to calculate the two variant QFIs for a pure 
state with different parameters. We adopt the standard 
notation where |1) = \\) and |0) = |t) correspond to the 
ground state and excited state, respectively. Consider an 
arbitrary single-qubit state 

|V^)=cos^|0) + e^^sin^|l), (28) 

of which the Bloch vector is denoted as u; = 
(sin ^ cos 0, sin ^ sine/), cos^)^, with and (j) referring to 
the polar and azimuth angles on the Bloch sphere. Here, 
the two parameters and (j) in Eq. (28) are assumed to 
be unitary encoded. 

Now we compute the QFIs of the single qubit state 
Eq. (28) in terms of the parameters: (amplitude pa- 
rameter) and (j) (phase parameter). With the help of Eq. 
(11), one can easily obtain To = 1 and = sin^ 0. For 
the amplitude parameter is constantly equal to 1, 

independent of both two parameters and (j). While 
about (j) depends on the parameter and reaches the 
maximum value 1, when = 7r/2, i.e., is a equal- 
weighted state = (|0) + e^^ |1)) /^/2 [10]. According 
to Eq. (25), we obtain Xq = 2 and = 2sin^ 6. Similar 
to Te^ Xq is also independent of the parameters and 
(j). Meanwhile, reaches the maximum value 2 at the 
point = 7^/2. In Sec. IIIB, we will assume that the 
qubit is initially in the equally weighted superposition of 
the two states |0) and |1), so both and reach their 
maximum values. 



III. QFI FOR A SINGLE QUBIT UNDER 
DECOHERENCE 

Decoherence occurs when a quantum system interacts 
with its environment, and it is unavoidable in almost all 
the realistic quantum systems. A quantum noisy dynam- 
ical process can be generally described by a map f , using 
the Kraus representation 

S{p) = Y,K^pKl (29) 

where are the Kraus operators satisfying 

Kj^ Ki_^ = U, which leads to the map £ being a 
CPT map [31, 32]. In the Bloch representation, such an 
non-unital trace-preserving process can be represented 
by an affine map on the generalized Bloch vector [34, 55]. 

For qubit system, Eq. (29) can be equivalently repre- 
sented as 

S{p) = ^l + ^{Au; + c)-&, (30) 

where A is a 3 x 3 real transformation matrix with el- 
ements defined as Aij = |Tr [a^f (a^)], and c G 
is the translation vector with elements given by q = 
lTY[ai£{l)]. From Eqs. (10) and (30), it indicates that 
under the decoherence, the Bloch vector u; of Eq. (10) 
is mapped as ^(cj) := A(jJ c. In parameter estima- 
tion, the unknown parameter is generally encoded on 
the probes through a unitary or non-unitary evolution 
[13]. In this paper, we do not consider the case of the 
non-unitary parametrization. A and c are assumed to 
be parameter-independent, i.e., the decoherence process 
will not introduce the parameter. 

With Eq. (11), the dynamic of the QFI based on the 
SLD operator for a single qubit under decoherence chan- 
nels can be generally expressed as 

^AHa.^(-)r + &^)^^^. (31) 
l-|£:(a;)| 

Similarly, according to Eq. (25), the dynamic of the vari- 
ant version of the QFI can be written as 

Xa = ^^£i^iL+e,(^) [5(a;) • d,S{u;)f , (32) 

Given an input state \iJj)^ the dynamics of the two QFIs 
under quantum channels are fully determined by the 
affine transformation matrix A and the translation vec- 
tor c. It is noted that Eqs. (31) and (32) are the general 
results that are applicable to all of those cases with dif- 
ferent parametrization processes. 

A. Dynamics of the QFIs under three decoherence 
channels 

Below, we will study the dynamics of the two vari- 
ant versions of the QFI under three paradigmatic types 
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of quantum channels [34, 37]: phase- damping channel 
(PDC), depolarizing channel (DPC), and generalized 
amplitude-damping channel (GADC) modeling a ther- 
mal bath at arbitrary temperature, which will be re- 
duced to the purely dissipative amplitude-damping chan- 
nel (ADC) when the environment temperature becomes 
zero. These channels are the prototype models of dissi- 
pation relevant in various experimental systems. 

Phase- damping channel: The PDC is a prototype 
model of dephasing or pure decoherence, i.e., loss of co- 
herence of a two-level state without any loss of system's 
energy. The PDC is described by the map 



= sp^p (poo|0)(0| + pii|l)(l|) , 



(33) 



and obviously the Kraus operators are given by 

KPDC = {V^1, Vp|0)(0|, (34) 

where p= 1 — s is the probability of the qubit exchanging 
a quantum with the bath at time t with s = exp (—7^/2), 
with 7 denoting the zero-temperature dissipation rate. 
For the PDC, there is no energy change and a loss of 
decoherence occurs with probability p. As a result of the 
action of the PDC, the Bloch sphere is compressed by a 
factor (1 — 2p) in the x?/-plane. 

According to Eq. (30), the transformation matrix A 
and the translation vector c of the PDC are given in 
Table I. These indicate that the Bloch vector components 
along the x- and ?/-axis shrink with probability 5, while 
the 2;-component remains invariant under the action of 
the PDC, and the Bloch vector uj is mapped as 

fpDC (<^) = {sUJ:r:,SUJy,UJzf • (35) 

Furthermore, with Eqs. (31) and (32), we obtain the ana- 
lytical results of the two variant versions of the QFI under 
the PDC in Table. I. One can find that J^^, Xq and are 
monotonic functions of t and are solely dependent on the 
parameter 0. Interestingly, Fe is constantly equal to 1 for 
any time, which implies that the QFI about the ampli- 
tude parameter 9 is robust under PDC. The significance 
of this is that one can avoid the impact of the PDC on 
the accuracy of the parameter estimation, by encoding 
the parameter on the amplitude of the input state. 
When = 7r/2, we have 

2 



To = 



1, J^d) 



1 

2 



2vT 



As is plotted in Figs. 1 (a) and (b), under the PDC, 
the QFIs associated to the phase parameter, and 
decreases monotonically with t and vanish only in the 
asymptotic limit t -> oo. Interestingly, Xg decrease from 
initial value 2 to final value 1, shown in Fig. 1 (b). It is 
indicated that the QFI Xq depend on the diagonal and 
off-diagonal elements of the density matrix. Due to the 
influence of the PDC, the off-diagonal ones vanish, and 
only the diagonal ones remain, which make Xq equal to 
1. 
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FIG. 1: (Color online) Plots of the two variant versions of the 
QFI versus t in terms of the parameter and under three 
quantum decoherence channels: PDC (a, b), DPC (c, d), and 
GADC (e, f) for a single qubit system. The qubit initially 
being prepared in the equally weighted state (Eq. (28) with 
^ = 7r/2) of which the QFIs take the maximum at t = 0. In 
(e) and (f), we plot the QFIs under the GADC for a zero- 
temperature n = and a finite-temperature n = 1, respec- 
tively. 



Depolarizing channel: The definition of the DPC is 
given via the map 



^DPc(p) = sp + p^. 



(36) 



and the corresponding Kraus operators are expressed as 



K 



DPC 



— 1, 



Vp Vp 



4 



(37) 



where s = 1—p and a = {ax, cTy, cFz) denotes the Pauli 
matrices. The channel is completely positive for < p < 
1. We see that for the DPC, the qubit is unchanged 
with probability s or it is depolarized to the completely 
mixed state 1/2 with probability p. It is seen that due to 
the action of the DPC, the radius of the Bloch sphere is 
reduced by a factor s, but its shape remains unchanged. 

The transformation matrix A and the translation vec- 
tor c for the DPC are given in Table I, and then the 
affine-mapped Bloch vector is obtained as 



^DPC {^) = {sUx, SUJy, SUz)^ 



(38) 



which shows that all components of the Bloch vector are 
shortened by a factor s. Moreover, we analytically derive 
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TABLE I: Analytical results for the time-evolutions of the two QFI quantities in terms of the parameters and (j) for the single 
qubit state in Eq. (28) under the quantum decoherence channels. The decoherence channels can geometrically be described 
as affine maps, and the transformation matrices A and the translation vectors c for each channel are listed. In order to simplify 

the expression of the results, we set X = sin^V^ 

= s[l + s+pcos (2(9)], and R2 =p^s'^ (a - /3 + cos (9)'' sin^ (9. 



Y = 2cos2 f ^s(l -s), Z = yjl-ssin^e- [p (a - /3) - s cos 6>]^ 



Quantum channel 



Phase-damping 
channel 
(PDC) 

Depolarizing 
channel 
(DPC) 

Amplitude Damping 
channel 
(ADC) 

Generalized Amplitude 
Damping channel 
(GADC) 
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s 




s / 
















s 


VI 
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-p(a-/3) 



sin^ e 



s sin 9 



^ + ^ ssin^^ 



3+4X + 2 cos^ 6'-cos(2 0) 
2(1+X)2 



2 - 2V1 - 



![3+4r+2 s sin^ ^+cos(2 6)] 
2(l+r)2 



1 + Z + ^2 (^^^ 



1+X 



2ssin2 6> 



1 + ^ 



the expressions of the two QFIs with respect to and 
under the DPC, given in Table 1. 

For 6> = 7r/2, those expressions are explicitly simplified 

as 



To 



Id) 



2vT 



The results show that, under the DPC, the two QFIs 
decrease by a factor of 5^ and s, respectively. The ex- 
pressions of the QFI T (X) for different parameters 6 
and (j) are the same. As are plotted in Figs. 1 (c) and 
(d), the QFIs decreases monotonically. 

Generalized Amplitude- damping channel: The 
GADC is given, in the Born-Markov approximation, via 
its Kraus representation as 



fGADc(p) = ^KipKl, 



(39) 



where the corresponding Kraus operators are 

Kgadc ={^/^ (|0) (0| + (1|) , ^/^|0) (1| , 



V^(^/i|0) (0| + |1) (1|), V^|l) (0|}, 



with 



and 



s, 



(40) 
(41) 

(42) 



where s and p are dependent on the mean number of 
excitations fi in the bath. In the zero-temperature limit, 
i.e., n = and s = s, Eq. (39) reduces to the purely 
dissipative ADC, and its Kraus operators are represented 
as 

Kadc = {\/5|0)(0| + |l)(l|, Vp|1)(0|}. (43) 

Similarly, in the Bloch representation, the GADC can 
be described as an affine map of which the transformation 
matrix A and the translation vector c being given Table 
I, and the Bloch vector uj is mapped as 

£gA-Dc{^) = (^V^UJx,V^^y,SUJz - p {a - P)^ . (44) 

When a = 1 and /3 = 0, Eq. (44) reduce to the ADC 
case. It indicates that the GADC squeezes the Bloch 
sphere into an ellipsoid and shifts it towards the north 
pole. The radius in the x?/-plane is reduced by a factor 
^/s, while in the 2;-direction it is reduced by a factor s. 
In the asymptotic limit t ^ oo, i.e., s = 0, p = 1, the 
Bloch vector becomes ^gadc (^) = (0, 0, — (a — /3))^, 
which also implies that, under the ADC (n = 0), the 
qubit finally stay in the ground state. Meanwhile, the 
analytical results of the two QFIs under the GADC (the 
finite temperature) and the ADC (the zero-temperature) 
are derived in Table I. 

When ^ = 7r/2, the dynamics of the QFI J^, under the 
ADC, in terms of and are the same, namely = 
Tcf) = s. As are shown in Figs. 1 (e) and (f), under 
the GADC and ADC, the QFIs with for the different 
parameters decrease monotonically with time. It is also 
shown that the QFIs for finite temperature decay more 
rapidly than that for zero-temperature. 
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B. Numerical calculation with hierarchy equation 

In this section, we focus on a simple dissipative model 
of a two-level system interacting with a zero-temperature 
bosonic reservoir to explicitly illustrate the behaviors of 
the two QFI quantities during time evolution [25, 38]. 
Here, we exactly examine this model by adopting the 
hierarchy equation method. The total Hamiltonian of 
the system and bath without performing the RWA is 



H 



1 



(45) 



with B = "^j^gkbk + h.c. Here, the first term is the 
free Hamiltonian of the qubit with transition frequency 
cjo, the second term denotes the environment part with 
the creation (annihilation) operators 6^. (bk) of the bath 
model with frequency cj/c, and the last term is the interac- 
tion Hamiltonian between the system and bath equipped 
with the coupling constant gk- In the zero temperature 
limit, the spectral density is generally represented by a 
Lorentzian [38, 39] 



1 



J(c^) = - 



A7 



{(jj - (jJo) + 7^ 



(46) 



where A reflects the system-bath coupling strength and 
7 is the spectral width of the coupling, which is related 
to the reservoir correlation time scale tb ~ 7~^. 

As is well known, this typical dissipative model is solv- 
able under the RWA which is effective in the weak cou- 
pling limit [38, 46]. With the RWA, the analytical time- 
evolution function of the system can be equivalently de- 
scribed as an ADC by redefining the Kraus operators 
Eq. (43) as 

Kadc = {h{t) |0)(0| + |1)(1|, Vi^-h' (t))|l)(0|} 



where h (t) is a crucial characteristic function as 



h{t) 



-7f/2 



(47) 



(48) 



with d = a/7^ — 4A. From Eq. (44), the afline-mapped 
bloch vector reads under this dissipative environment 

^ADC {U^) = {h {t) UJa:, h {t) UJy, 1 - {t)f . (49) 

As shown in Table I, then, one can easily recover the re- 
sult of the dynamical QFI in terms of given in Ref. [25] , 



^0 = h\t) . 



(50) 



With the hierarchy equation method, the exact dy- 
namic of the system is derived as the following equation 
in the interaction picture [46] 

p^^) (t) = r exp [- j^dt2 dhV {t2T [C^ {t2 - h) 
V (ti)^ + iC' (t2 - h) V (ti)° ] }ps(o), (51) 



where T is the chronological time-ordering operator and, 
to simplify the description, we introduce two super- 
operators A^B = [A,B] and A°B = {A,B}. Assume 
that the initial state is taken as p (0) = ps (0) pB with 
Pb being in the vacuum state 0/c |0/c). Also, (^2 — h) 
and (^2 — ^i) respectively correspond to the real and 
imaginary part of the bath time-correlation function 
which is defined as 

C {t2 - h) = {B {t2) B (ti))B = Tr [B (^2) B {h) pb] , 

(52) 

where B {t) = Y.^gkhue'''^^^ + h.c. Then the time- 
correlation function Eq. (52) becomes the exponential 
form 



C{t2 - ti) = A exp [- (7-hia;o) Ih - ti\ 



(53) 



With Eqs. (51) and (53), we further obtain the set of 
hierarchical equations of the qubit as [46] 

d ^ 

k=l 

- 4i2^k\v-+i-lfV°]gH.g,{t), (54) 



k=l 

where the subscript n = (ni, 77,2) is a two-dimensional 
index with 71^2) ^ 0? and ps {t) = ^(0,0) (^)- The vec- 
tors are ei = (1, 0), 62 = (0, 1), and P = (z/i, 1^2) = 
(7 — iujQ, 7 -h iuJo). We emphasize that Qfi (t), with n 7^ 
(0, 0), are auxiliary operators introduced only for the 
sake of computing, they are not density matrices, and 
are all set to be zero at t = 0. Through solving the above 
hierarchy equations, the dynamics of the system can be 
exactly determined without making the RWA. 

Figure 2 displays the dynamics of the QFI of Eq. (2) 
and the variant QFI of Eq. (16) versus time for differ- 
ent parameters initially encoded in the qubit state. The 
spectral width of the coupling is set by 7 = 0.2co'o- In 
the Markovian regime (inset plots with A = O.OI7), the 
QFIs J^0((f)) and ^6'(0) ^ire monotonically go to zero. Both 
the numerical and the analytical results are consistent 
in the weak coupling limit. While, in the non-Makovian 
regime (7 = O.IA), due to the strong couphng with the 
reservoir, the dynamics of these two information quanti- 
ties exhibit the oscillations and revivals over time. The 
times when vanish completely coincide with those 

for J-6>(0). The surprising aspect here is that the quan- 
tities of ^0((f)) for the two parameters and (j) decrease 
quickly to zero from the initial value of 2, and then 
later on revive with almost the same amount as J-'e{(f>) 
do. Since the hierarchical equations we obtained exactly 
depict the dynamics of the system, the oscillating devi- 
ations between the numerical and the analytical results 
can be interpreted as the contribution of the counter- 
rotating- wave terms which have been omitted by making 
the RWA [46] . Such deviations behave evidently for both 
two information quantities in terms of ^. As shown in 
Ref. [25], they addressed these oscillations and revivals 
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FIG. 2: (Color online) Two QFIs Eq. (2) (a, b) and Eq. (16) 
(c, d) versus time for different parameters initially embedded 
in the qubit under the dissipative process. The solid lines 
display the numerical results by using the hierarchy equation 
method. These are compared to the analytical results with 
RWA (dashed lines), all plotted in the non-Markovian regime 
(A = O.I7). The insets show in the Markovian regime (A = 
O.OI7) 



in Figs. 2 (a) and (b) by introducing the definition of 
the QFI flow. Analogously, we also introduce a variant 
QFI flow (or skew information flow) as the change rate 
of information quantity X, i.e., a := dXx/dt. Also, a > 
indicates that there is a information flow from the en- 
vironment to the system, corresponding to the revivals, 
and cr < denotes that the flow from the system to 
the environment, accounting for the decays, as shown in 
Figs. 2 (c) and (d). 



where Tr (p^) is the purity. The purity equal to one, 
corresponding to pure states, and l/d, corresponding to 
mixed states. With Eq. (55), we get 



1 1 , , 



Tr (p^) 

by using following relation 

Tr [(a ■ 17) (6 • 17)] = 2 a ■ 6. 



(57) 



(58) 



Thus, from Eqs. (56) and (57), we obtain the length of 
the generalized Bloch vector satisfying 



< a; < ^2 {d - 1) /d, 



(59) 



where the first equality holds for the maximal mixed state 
and the second for pure states. 

Owning to the completeness of the generators of Lie 
algebra, any Hermitian matrix can be described by the 
common set of generators. Based on this representation, 
we give the new expressions of the two QFIs in terms of 
the Bloch vector. Firstly, we consider the QFI in Eq. (2). 
Tx can be represented as 



{dxu;)' M-^dxu;, uj < ^2 {d - 1) /d, 



\dxuj\ 



UJ = V'2 {d - 1) /d, 



where is a real symmetry matrix defined as 
2 
d 



(60) 



(61) 



The superscript T in the above equations denotes the 
transpose operation and denotes the matrix inverse 
of A4. 11^2_i is the identity matrix of dimension d^ — 1 
and G is a (d^ — l) x (d^ — l) real symmetric matrix 
whose ij-element is 



IV. GENERALIZATION 

A. Representation of the QFI in terms of the 
generalized Bloch vector for a qudit system 

We have examined the QFIs under decoherence for a 
single qubit in the Bloch representation. Now, we con- 
sider a more general case, i.e., the qudit (a d-dimensional 
quantum system). A general qudit state can be written 
in the Bloch representation as [56, 57] 



\G] 



P 



1 

2 



(55) 



where 1^^ is a d x (i identity operator, i) = {rji}^^^ are 
the generators of the Lie algebra 5U (d) (see appendix 
A), and u; G ~^ denotes the generalized Bloch vector 
whose ith element is Tr(p7)^). 

Since we have the following relation 



1/d < Tr (p^) < 1, 



(56) 



k=l 



(62) 



where gijk is the completely symmetric tensor defined 
in Eq. (A. 7). Hence, A4 is also real symmetric matrix. 
Since A4 may have some zero eigenvalues, the inverse is 
defined on the support of M. [55]. Here, the first line of 
Eq. (60) only applies to mixed states, and the detailed 
derivation can be found in appendix B. 

For pure states, Tx is generally expressed as the norm 
of the derivative of the Bloch vector shown in Eq. (60). 
It can be easily derived by using Eq. (58) and following 
the same procedure used in deriving the second line of 
Eq. (11). 

From Eq. (61), we can see that the matrix M is de- 
pendent of UJ. Then we can conclude that J^x given by 
Eq. (60) only depends on the two elements: the Bloch 
vector of the density matrix and the derivative of it. 
Moreover, it deserves to emphasize that Eq. (60) is also 
valid for arbitrary quantum state, since a density ma- 
trix always can be expressed as the form of Eq. (55) by 
expanding over the generators of the Lie algebra. 
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When d = 2, we have G = 0, due to gijk = 0. Then, 
the real symmetric matrix Ai reduces to 



M = Is - (jJijS^ 
The inverse of M is verified as 



M- 



13 + 



1 



(63) 



(64) 



Substituting the above equation into Eq. (60) finaUy re- 
covers the first line of Eq. (11). 

We next take account of the QFI in Eq. (16). We firstly 
expand 



Then, X\ can be represented as follows 



(65) 



{dxyf + l\dxx\' , uj<^2{d-l) /d, 



oj = y/2 {d - 1) /d, 



(66) 

where y and x are completely determined by the follow- 
ing (P quadratic nonlinear equations: 



2,^1 |2 

y 



1 



9 -L — ^ 

^ y ~r / ^ Qijk Xj — ^ , 



(67) 
(68) 



for /c = 1, 2, — 1. The first line of the above equation 
is applicable for mixed states, and the detailed derivation 
is given in appendix C. 

For pure states, Xx is generally expressed as the norm 
of the derivative of the Bloch vector up to a factor of 2 
shown in Eq. (66). It can be obtained by using Eq. (58) 
and following the same procedure used in deriving the 
second line of Eq. (25). 

For the case of d = 2, equations (67) and (68) reduce 

to 



y 



^\x\ 
Ayx 



1 

2' 



(69) 
(70) 



with Qijk = 0. By solving the above equations (see ap- 
pendix C), one can obtain 



y 

X 



2 



1^1 + ^ 



(71) 
(72) 



Inserting the above solutions into the first line of Eq. (66) 
and making some simplification recovers the first line 
given in Eq. (25). 



B. QFI for an A/^-qubit system in noisy environment 

Below, we study the dynamic of the QFI for an N- 
qubit system in noisy environment. In the Bloch repre- 
sentation, the two QFIs for qudit are described in terms 
of the generalized Bloch vector a;, as is shown in Eqs. (60) 
and (66) respectively. The Bloch vector is assumed to be 
the function of an unknown parameter A on the system. 
Meanwhile, we emphasize that Eqs. (60) and (66) are 
also applicable for the multi-qubit system with exchange 
symmetry. Since, a collection of TV qubits is represented 
by the collective operators [58] 



N 



,{a = x, y, z) , 



(73) 



where (Jia denotes the pauli matrix of the ith qubit. 
Such an A^-qubit ensemble with total angular momen- 
tum j = N/2 can be approximately viewed as a qudit 
system, when it has the symmetry under the exchange of 
two qubits. The collective basis of this system is m)} 
for m = 0, dzl, dz2, itj, which is so-called Dicke state, 
Jz \ j^m) = m|j, m). Hence, Eqs. (60) and (66) are also 
valid for those multi-qubit systems. 

Assume that the dimension of the decohered state £{p) 
is the same of that of p. Similar to Eq. (30), an 7V-qubit 
system with exchange symmetry described by Eq. (55) 
under decoherence can be expressed as [34, 55] 

e{p) = -^td + \{Au: + c).'n, (74) 

with d = N Here, A is a matrix of dimension d^ — 1 
with elements 



Ai^ = -TY%e{fi,)], (75) 
and c is a c/^ — 1 dimensional vector with elements 



IV [f)i S{ Id)] ■ 



(76) 



Equation (74) illustrates that a Markovian quantum dy- 
namic can be geometrically described as an affine trans- 
formation, i.e., 



£ : uj A u; -\- c. 



(77) 



As we mentioned at the beginning of Sec. HI, in the pa- 
rameter estimation, the unknown parameter is generally 
imprinted into the probes through unitary or non-unitary 
evolution [13]. It is noted that equation (60) is applicable 
to the cases with different the parametrization processes, 
by replacing the Bloch vector cj with the affine-mapped 
Bloch vector £{otj). 

To be more specific, we consider an experimentally re- 
alizable Ramsey interferometry [59], to estimate a phys- 
ical parameter in a Bose-Einstein condensate (BEC) of 
N two-level atoms interacting with a common thermal 
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FIG. 3: (Color online) The Schematic representation of Ram- 
sey interferometry in presence of collisional dephasing. 



reservoir. This experiment we model is shown schemat- 
ically in Fig. 3. Those two- level atoms may be con- 
sidered as qubits or probes, which are prepared in an 
A/'-qubit Greenberger-Horne-Zeilinger (GHZ) state (or 
Schrodinger-cat state) 



GHZ/ 



i=(io>«-+iir). 



(78) 



Considering all the qubits being initially prepared in |0), 
then a GHZ state is generated by putting a Hadamard 
gate acts on the first qubit, followed by a sequence of 
controlled-NOT gates linking the first one with each of 
the remaining ones [10, 12, 59]. It is well known that 
such a state saturates the ultimate Heisenberg limit (HL) 
1/N on precision of measurement. This result gives a 
quadratic improvement over the standard quantum limit 
(SQL), which is achievable with product states [10, 59]. 

As is shown in Fig. 3, the parameter cp of interest is 
unitary imprinted on the state of the probe qubits. The 
unitary operator is given by 



u: 



tot 



exp 



(79) 



where Jz is the z component of the total angular momen- 
tum for all qubits [59]. After the unitary evolution, then 
the state of the probes becomes 



V^GHZ 



1 

1 

71 



(|0)^ 



N N 
2"''2 



il 11 \ -Lpi^V 



?,-r)).(so) 



up to a global phase factor e ^^'^1^ , The second equal- 
ity above is valid in the standard representation of the 



generator J^. According to the appendix D, the Bloch 
vector of the state of Eq. (80) reads 



fcos (7V(/9) , 0, 0, sin {Nip) , 0, 0, 



id^-d)/2 
1 



{d?-d)i2 



^/2m {m -\- 1) 
2-d_\T 



V2 {d-i){d-2y 



y/2d{d-l)^ 



(81) 



for m = 1, (i — 2. For the sake of clarity, we take N = 2 
(i.e., d = 3) for example. The Bloch vector for N = 2 in 
Eq. (81) becomes 

1 It 

u;^ = {cos{2(p) , 0, 0, sm{2(p) , 0, 0, -, -^-^) ^ 

of dimension 8. As is shown in Eq. (81), u;^p only con- 
tains two (/^-dependent elements. From Eq. (60), we find 
that that those (/^-independent elements of uj^ will not 
contribute to computation of J^^. 

In the realistic experiment, decoherence always exists. 
As is shown in Fig. 3, we consider the effect of the col- 
lisional dephasing on this measurement protocol, which 
is induced by the interaction between the qubits and the 
common thermal reservoir [60, 61]. The master equation 
of the system in the Lindblad form can be described as 
[61, 62] 

p{t)=Cp^j (2Jzp it) Jz-p it) Jl - JIp it)) , (82) 

where C denotes the Lindblad superoperator, 7 denotes 
the dephasing rate, and p is the reduced density operator 
of the system in the interaction picture. For the single 
qubit case, Eq. (82) reduces to 



p{t) = Cp=l[cjzp(t) cJz-p{t)], 



(83) 



which corresponds to a single-qubit dephasing channel 
(i.e., DPC). Here, we use the interaction representation 
which does not affect the result of the calculation, since 
the QFI remains invariant under the unitary evolution 
being independent of the parameter (p. From Eq. (82), 
the time evolution of the density matrix elements is given 
as follows 

Pm,n it) = {j, m\ p (t) \j, n) = pm,n (0) e"'™-")'^*. (84) 

In the Bloch representation, the corresponding affine 
transformation matrix for the collisional dephasing in 
Eq. (82) can be obtained as 
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Diag(^e 



d-l 

P N^-ft ^-{N-lf-ft ^-{N-lf-ft -iN-2f-ft —ft —ft 



((i2_d)/2 



d-l 



N^7t^^-{N-ir^t ^-{N-iy^t ^-{N-2Y-ft ^-A-ft ^-^t 



,e 

(d2-d)/2 



, . . . , 6 ,6 



e"^^ 11 11 

,...,0 , X, X, -LJ? 



(85) 



d-l 



and c = (see appendix D). Apparently, both transfor- 
mation matrix A and translation vector c are indepen- 
dent of ip. After the dephasing process, the Bloch vector 
of Eq. (81) is affine-mapped as 

8 (cj^) = (e-^'^^ cos [N^) , 0, 0, e'^'^' sin [N^) , 



0.09 



a/2 (d - 1) (d - 2) ' V2rf ((i- 1) 



(86) 



which is given by inserting A of Eq. (85) and c = into 
Eq. (77). 

The dynamic of the QFI associated to the parameter 
may be evaluated as, from Eqs. (60) and (86), 



-2N^-ft 



(87) 



which shows that the quantity of the QFI is monoton- 
ically decreased by a factor of e~^^ The ultimate 
accuracy of the estimation reads 

= J^^^,, (88) 

by inserting Eq. (87) into Eq. (6). It indicates that the 
collisional dephasing dramatically deteriorate the sensi- 
tivity of the phase (p estimation. Our result is consistent 
with Ref. [63]. 

To clearly see the effect of the decoherence process, we 
define a time scale tc which is the time over which the 
QFI reduce from the HL (i.e., J^^ = N'^) to the SQL 
(i.e., = N) [6]. For the collisional dephasing, the 
characteristic time reads 

log .X 
= (89) 

by setting 7 = 1. After above, the GHZ state cannot be 
used to perform over shot-noise estimation. As is plotted 
in Fig. 4, the characteristic time decreases exponentially 
as N increases. It illustrates that the advantage of a GHZ 
state deteriorates in the case of collisional dephasing. 



V. DISCUSSION AND CONCLUSION 

We have discussed the time evolution of the two variant 
versions of the QFIs in the presence of quantum noises. 



0.06 



0.03 




FIG. 4: (Color online) Plot of the characteristic time tc in 
Eq. (89) (blue dash-circle curve). 



With the help of the Bloch representation, we derived 
the explicit formula of the two information quantities for 
a single qubit system. The analytical expressions of the 
dynamics of the two QFIs under three typical quantum 
decoherence channels were obtained. Both information 
quantities in those channels are decreased monotonically 
with time, only except the case that the J^q remained 
invariant under PDC. It manifested that the QFI defined 
in Eq. (2) about the amplitude parameter 9 is robust for 
the PDC. 

We also considered a simple dissipative model of a sin- 
gle qubit coupling with a bosonic reservoir at zero tem- 
perature. By applying the hierarchy equation method, 
we exactly calculated the dynamical QFIs during time 
evolution. We found that the numerical results qualita- 
tively coincided with the analytical ones using the RWA. 
The deviation between them is accounted as the contri- 
bution of the counter-rotating- wave terms. In the weak 
coupling regime, we observed that two QFIs about dif- 
ferent parameters were monotonically decreased. When 
the strength of the coupling become more stronger, the 
behavior of the non-Markovian could be observed from 
the QFIs perspective. 

Finally, we generalized the results to the qudit system, 
expressing the two QFIs in terms of the generalized Bloch 
vector. Those expressions were valid for the A/'-qubit sys- 
tem with symmetry exchange. We also considerd the QFI 
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in the presence of the colhsional dephasing with the ini- 
tial state being prepared as a GHZ state. The affine 
matrices for this dephasing process was derived and the 
Bloch vector of the GHZ state also was expressed as the 
Bloch vector. We found the dynamical QFI exponentially 
decrease by a factor of e~^^ 



Employing completeness relation of generators of 5u((i), 
one can expand an arbitrary dimensional Hermitian 
matrix X as 



X = - Tr (X) Id + ^ Xk fjk, Xk e 



(A.8) 



k=l 
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Appendix A: Generators of the Lie algebra su (d) 

In group theory, the generators of the Lie algebra have 
the following properties: i) Hermitian: 



ii) traceless: 



Trr), =0; 



(A.l) 



(A.2) 



iii) orthogonality and normalization with respect to the 
trace metric relation: 



-Ti{fiif]j) =5ij. 



{A.3) 



Moreover, they also satisfy two relations as follows, 
characterized by the structure constants fijk and gijk, 



[fii, fjj] = 2i^ fijkVk, 



(A.4) 



{Vi, %} = -^^ij^d-^-'^^ gijkVk, (A.5) 



where 1^ is the unit matrix of dimension and fijk 
{dijk) denotes the completely antisymmetric (symmetric) 
tensor. The structure constants are determined by 



fijk = ^Tr([i7i,%]i7fe), 

9ijk = ^Tr({7}i,%}7}/e). 



(A.6) 
(A.7) 



with Xk = Tr (Xry/e). 

Below, we systematically construct the generators fj = 

{Vj}j=i of ^he Lie algebra 5U (d) which are given as fol- 
lows [56, 57]. On a d-dimensional Hilbert space H G 
spanned by an orthonormal set of states {| Ta)}^^-^^, we 
first construct two sets of block off-diagonal Hermitian 
traceless matrices, 



\m) {n\ -h |n) (m| , (A.9) 
-i{\m){n\-\n){m\), (A.IO) 



for 1 < m < n < (i. The above two sets contain the same 
number of elements as (^d^ — d^ /2. Then, we construct 
d — 1 real diagonal traceless matrices 



Vk 



k{k- 



k 

-[^ |m) {m\-k |A: + 1)(A: + 1| 

^ m=l 



(A.11) 

for {1 < k < d — 1). Hitherto, the generators for the Lie 
algebra 5u{d) are obtained by the new set [56, 57], 



{Vi}f=i^ = { 



^m,n^ ^m,n^ ^k 



}■ 



(A.12) 



The generators for (i = 2 coincide with the Pauli matri- 
ces with the structure constants fijk being Levi-Civita 
symbol eijk and gijk = 0. The generators of su (3) are 
described by the Gell-Mann matrices given by [57] 



m 



V3 



V5 



V7 




V6 




V8 



1 
75 I 1 I . (A.13) 
* -2 



Appendix B: Derivation of equation (60) 

Now, we provide the detailed derivation of the first line 
of Eq. (60) . In order to express in terms of the Bloch 
vector. We should firstly write all the Hermitian opera- 
tor in the same representation, expanding them over the 
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common generators of the Lie algebra 5u{d) [55]. The 
density matrix p is described by Eq. (55), and the SLD 
operator is supposed to be expanded as 



L = aid ^b rj, 



(B.l) 



where a and b are respectively real number and real vec- 
tor to be determined. 

Suppose p contains some unknown parameter A, i.e., 
the Bloch vector u; is A-dependent. Then we have 



1 T 
9xP = 2 (^Au;) i). 



(B.2) 



By substituting Eqs. (B.l) and (B.2) into Eq. (2), we find 
that J^x can be rewritten as 



Tx = {dxu;fb. 



(B.3) 



Here, we use Eq. (58) and the traceless property of gen- 
erators. Equation (B.3) shows that when b is derived, 
then the problem is solved. 

To determine 6, we should use Eq. (3). Furthermore, 
by using Eq. (A. 5), we have following relation 

{a • r), 6 • r)} = ^ aibj {fji, %} 



- a^b Irf + 2 ^ QijkCiihj r}/c, (B.4) 

ijk 



With the help of Eqs. (55), (B.l), and (B.4), the right- 
hand side (RHS) of Eq. (3) reads 

\{P.L} = i({i,aIL + 6.r7} + 

{-a? -17, a]l} + {-a;-i7, 6-17}) 

= i (a + u;^b) 1 + 

k ij 

The left-hand side (LHS) of Eq. (3) is given by Eq. (B.2). 
By comparing the terms on both side of Eq. (3), one 
obtains 

a^uj^b = 0, (B.5) 
{dxu:)^ 7) = ^ (^-bk + aujk + ^gijkOJihjjr]k 

k ij 

2 

= -^^^r) + auj^T] + ^ gijk^ihjfjk 

ijk 



jk 

b^rj^acj^Tj^b^Grj. 



Here, the matrix element of G in Eq. (B.6) is given by 
Gjk = Y.i9ijk^i satisfying Gjk = G/cj. With Eq. (B.5), 
we have 



a = —uS^b = —b^uj. 



(B.7) 



By insetting the above equation into Eq. (B.6), we obtain 



{dxujf r) = 6^ ( 3 - cjcjT + G ) r). 



by setting 



b^Mfi, 



M = -- uu;^ + G. 
a 



(B.8) 



Hence, equation (B.8) directly gives the following equa- 
tion 



(B.9) 



Suppose that we have the inverse matrix . After 
Eq. (B.ll), we will make further discussions about . 
We have 



h = M-^ {dxu;) . 



(B.IO) 



from Eq. (B.9). Finally, inserting Eq. (B.IO) into 
Eq. (B.3) yields 



(B.ll) 



To calculate J^a, we need to find the inverse of M. 
Generally, it may or may not exist, i.e., A4 may have 
some zero eigenvalues. In this case, we define A4~^ on the 
support of i.e., supp(Al), which is defined as a space 
spanned by those eigenvectors with nonzero eigenvalues. 
It is reasonable to do so. By inserting Eq. (B.9) into 
Eq. (B.3), we rewrite Tx 



= h^Mb. 

Suppose M. has spectral decomposition 

i^l 



(B.12) 



(B.13) 



where Vi denotes eigenvector with eigenvalue rrii. We 
assume rrii for i = 1, n and = for i = n -\- 
1, ...,^2 - 1. With Eq. (B.13), equation (B.12) reads 



= y^^mib^VjvJb 



(B.14) 



(B.6) 

It is shown that J^x is defined on supp(Al). 



14 



Appendix C: Derivation of equation (66) 

In this appendix, we give the detailed derivation of X\ 
in terms of the Bloch vector for mixed states given by 
Eq. (66). We first expand the Hermitian matrix as 



(C.l) 



where y = Tr (y^) /d and x is an unknown — 1 dimen- 
sional real vector. With Eq. (C.l), I\ in Eq. (16) can be 
rewritten as 



{dxyf + I \dxx\^ 



(C.2) 



by using the properties of the generators of the Lie alge- 
bra given in App. A. Furthermore, we have the following 
equation 



(C.3) 



In the Bloch representation, the left-hand side of the 
equation above is given by Eq. (55), and the right-hand 
side reads 

{Vpf ^ {^^ ^ \ 1^1^) ^d^-i + + QijkXiXj) fjk. 

By comparing the terms on both sides, one finds that 
y and x are completely determined by the following dP 
quadratic nonlinear equations: 



2 2 I i2 I 



^ ^ Qijk Xi Xj -\- 1y X]^ 



2 ' 



(C.4) 
(C.5) 



for k = 1,2,..., dp - 1. 

Below, we give the detailed derivation for Eqs. (71) 
and (72). Multiplying the LHS of Eq. (70) hy Ayx^ and 
the RHS by a;^, we obtain 



-\a 2| |2 I |2 

my \x\ = |a;| . 



By replacing \x\^ in equation above by 



1 



\x\ 



(C.6) 



(C.7) 



given by Eq. (69), we get the following equation 



16/-8^^ + |u;|^ = 0. 



(C.8) 



Solving the above equation, we obtain two solutions 



(C.9) 
(C.IO) 



We suppose that p has spectral decomposition 

P = ^Qi = ^1 + (1 - ^i) l'02) {'02| . 

(C.ll) 



Then we can derive the following inequality 



(C.12) 



Then the solution y'^ in Eq. (C.IO) can be neglected. 
Due to y being nonnegative real number, we finally 
get Eq. (71) from Eq. (C.9). Furthermore, substituting 
Eq. (71) into Eq. (70) directly yields Eq. (72). 



Appendix D: Derivation of the afRne transformation 
matrix A and c for the collisional dephasing 

We firstly verify the elements of the Bloch vector can 
be directly written out from the elements of the density 
matrix. A d-dimensional density matrix is described by 
the generalized Bloch vector in Eq. (55) in the Bloch 
representation. The elements of the generalized Bloch 
vector is defined by 



(D.l) 



In the appendix A, we systematically construct the gener- 
ators of Lie algebra su {d) in Eq. (A. 12), which is defined 
by three sets given by Eqs. (A.9), (A.IO), and (A.ll). 
Then a; also can be divided into three parts as 



r W^-i r 
{^i}i=i ={ 



^ A 



n5 



(D.2) 



We find that the elements of the Bloch vector are directly 
given by the elements of the density matrix 



Tr (pi™,„) = -2Q{pm,n), 



(D.3) 
(D.4) 



Tr 



k 

^ ^ m=l 



kpk+i,k+i -(0.5) 



The master equation of the collisional dephasing model 
is described by Eq. (82), and the time evolution of the 
density matrix elements is given by Eq. (84). With 
Eqs. (84), (D.3), (D.4), and (D.3), the elements of the 
affine-mapped Bloch vector is given by 

<n(*) = e-(— )'^*2KK,„(0)] 

= e-(— ")'^*c.^,„(0), (D.6) 
_g-(™-„)V2c>K,„ (0)] 



-{m—n) 7t A 



,n(0). 



(D.7) 
{D.8) 
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Equation (D.8) illustrates that u;^ remains unchanged changed by dephasing. Writing in the form of affine map 
under the collsional dephasing, since lv^ depend on the in Eq. (74), we finally obtain A in Eq. (85) and c = 
diagonal elements of the density matrix which left un- from Eqs. (D.6), (D.7), and (D.8). 
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